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Abstract. A method is presented here for investigating variations in the upper end 
of the stellar Initial Mass Function (IMF) by probing the production rate of ionizing 
photons in unresolved, compact star clusters with ages<10 Myr and covering a range 
of masses. We test this method on the young cluster population in the nearby galaxy 
M51a, for which multi-wavelength observations from the Hubble Space Telescope are 
available. Our results indicate that the proposed method can probe the upper end of the 
IMF in galaxies located out to at least 10 Mpc, i.e., a factor 200 further away than 
possible by counting individual stars in young compact clusters. Our results for this 
galaxy show no obvious dependence of the upper mass end of the IMF on the mass of 
the star cluster, down to ~1000 M , although more extensive analyses involving lower 
mass clusters and other galaxies are needed to confirm this conclusion. 



1. Introduction 

Recent analyses of optical and UV data of nearby galaxies have clearly indicated that, as 
some global galactic parameter decreases (luminosity, SFR, SFR/area, etc.), so does the 
galaxy-integrated ionizing photon rate, as traced by the hydrogen recombination line 
Ha, per unit UV or o ptical luminosity dHoversten & G lazebrook 2008]; iMeurer et al. 



Ha, p er unit UV or optical luminosity ( Moverste n & utazebrooKl zUUo; Meure r et at. 
1200 9|; iLee et al.1l2009l) . Presentations from these authors, containing also updates over 



the cited papers, are included in these Proceedings, together with similar results ob- 
tained by other investigators. 

In all cases, corrections for dust attenuation have been included and/or discussed, 
and uncertainties deriving from such corrections are unlikely to produce the observed 
trend over the full luminosity range analyzed (factor *10 4 in B-band luminosity, see 
ILee et al.1l2009 ). At least three possible, non-mutually-exclusive, scenarios have been 



proposed to account for the observed trend: 
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1. it is an effect of stochastic sampling and/or a particular period in the star for- 
mation history of the ionizing photon rate-deficient (N, on -deficient) galaxies. 
Specifically, these galaxies could be in a post-burst phase, and the apparent 
higher frequency of such phase in low-luminosity systems could be due to the 
sporadic nature of star formation in these systems; 

2. it is an effect of ionizing photon leakage from the star-forming regions of these 
low-mass and low-density systems; 

3. it is an effect of a systematic behavior of the high-end of the stellar IMF linked 
to the variation of some global galactic parameter. 

The first scenario has been discussed by all authors cited so far, and considered 
unlikely to be the only or primary driver of the observed trend, given that it would re- 
quire synchronization among the star formation histories of the low-luminosity galax- 
ies (see, also, the presentation by E. Hoversten in these Proceedings). These results 
notwithstanding, it may still be worth investigating the potential effects of sporadic 
star formation with increasingly sophisticated prescriptions for the star formation his- 
tories of galaxies. Accurate star formation histories are only now becomin g available 



for a significant number of low-mass and low-luminosity galaxies (e.g., IWeisz et al. 



2008). Furthermore, low-luminosity galaxies show over an order of magnitude large r 
variation in their birth-rate parameter than higher luminosity galaxies ( Lee et al1l2007 ). 

which may possibly impact observables. 

The second scenario is discussed in, e.g. jHunter et all ( |2010|) . Although a poten- 
tially attractive interpretation, leakage of photons from galaxies is difficult to measure. 
A few observations performed in higher density systems than those under considera- 
tion here place an uppe r limit of a few percent to the leakage of ion i zing photons from 
low-redshift galaxies (iHeckman et al.l l 200ll ; fDeharve ng et al.ll200ll : lHayes et ai"1l2007 



Bland-Hawthorn & Malonevl ll999ir Similar limits of about 6%, and no higher than 



20%, have been recently obtained for Lyman Break Galaxies at redshift «4 (IVanzella et al 
201dh . Ionizing photons leakage from low-luminosity, low-mass galaxies is clearly im- 
portant to establish, as it could have major implications for the re-ionization of the Uni- 
verse at high-redshift. However, at low density, leakage of ionizing photons throughout 
a galaxy disk or into a diffuse halo would be very difficult to detect at the resulting low 
emission measures, even if the photons did not l eave the galaxy entirely. 

The third scenario has bee n suggested by iPflamm-Altenburg et all d2009h. and 
builds o n the model propose d bvlWeidner. Kroupa & Larse nl d2004l).T Weidner & Kroupa 



(120051) . IWeidner & Kroupal <|200eJ), and Weidner, Kroupa & Bonnell (2010|), accord 
ing to which: (1) there is a correlation between a galaxy's S FR and the maximum 
mass of the cluster that can form in that galaxy dLarsenll2002l) : (2) there is a correla- 
tion between a cluster ma ss and the most massive possible star the cluster can form 
(Weid ner & Kroupal 120061) . The combination of these two conditions produces an in- 
tegrated galaxial stellar IMG (IGIMF) that reproduces the observed decrease in the 
ionizing photon rate per unit UV luminosity as a function of decreasing galactic SFR, 
if the assumption that most of the stars form in stellar clusters holds. 

1.1. Goal of This Presentation 

This presentation, based on the results in ICalzetti et all ( 2010h . is concerned with the 
second of the IGIMF's requirements, i.e., that a correlation might be present between a 
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cluster mass and the most massive star it can form. This is another way to say that the 
form ation of stars in a cluster is not an event solely driven by stochastic sampling of the 
IMF dElmegreerfl2006T) . and 100 clusters each with mass=10 2 M G are not equivalent to 
one 10 4 M Q cluster, as the former will systematically lack the high-mass stars contained 
in the latter. 

Measuring the number of massive stars in a stellar cluster is a problem that in- 
creases in difficulty with increasing distance. Signific ant blending will occur for mas- 
sive s tars due to crowding in the centers of star clusters (Mafz-Apellamz 2008; lAscenso et aL 



12009b . even with the angular resolution of the Hubble Space Telescope, for distances 
larger than those of the Magellanic Clouds (~55 kpc). All the N; „-deficient galaxies 
are farther away than the Magellanic Clouds. Furthermore, the labor-intensive ap- 
proach of counting individual stars in single-age, young stellar populations (typically 
<3 Myr old star clusters), which is the standard approach for deriving stellar IMFs, 
cannot build statistically large samples of such measurements. 

The approach to constraining the upper end of the IMF presented here exploits the 
integrated properties of star clusters. Each star cluster is treated as a single unit, and 
many such units are averaged together in order to investigate variations on the upper- 
IMF; this is accomplished by measuring the dependence of the ionizing photon rate, 
normalized to the cluster mass, on cluster mass. This requires t hat accurate ages and 
masses be determined for each cluster. This method extends the Corbelli et all (2009) 
approach (the cluster birthline, see also Corbelli 's contribution to these Proceedings) 
by normalizing the ionizing photon rate to the age-independent cluster mass, instead of 
the age-dependent bolometric luminosity. 



2. Framework 

Measuring the ionizing photon rate, N !0 „, in a young, massive star cluster is equivalent 
to measuring the number of massive stars that cluster contains at that age. The younger 
the age, the closer the measured number is to the total ionization rate of massive stars at 
birth (including the most massive stars formed), or the integral of the upper-IMF over 
mass. Star clusters also represent ideal sites for such measurements, because massive 
star populations tend to be coeval or close to coeval. If a trend is predicted for the 
number of massive stars or the most massive star formed as a function of the cluster's 
mass, then the presence/absence of such trend can be established by investigating the 
behavior of N,- on /M c / as a function of the cluster's mass M f /. The left panel of Figure Q] 
shows the expected trends for a 'universal' IMF, where the most massive star formed is 
independent of the cluster's mass, and for an I MF where the most massi ve star formed 
depends on the cluster's mass (the model of IWeidner & Kroupal 12006): the ionizing 



photon rate is here replaced with the extinction-corrected Ha luminosity. The Figure 
also shows the expected range for different values of the cluster's metallicity and for 
different ranges of the clusters' ages. 

This simple framework is fraught with a number of uncertainties, all of which need 
to be carefully evaluated. Even in the event that the IMF is 'universal', stochastic sam- 
pling of the IMF is one of such uncertainties, which starts to be important as the mass of 
the cluster decreases. Below masseswa fewxlO 4 M Q , the IMF cannot be fully sampled, 
and since fractional stars cannot be created, the stars formed in the cluster will be ran- 
domly sampled from the IMF. This induces a scatter on the measured parameters. For 
instance, the scatter on the mass determination increases from 6% for a 10 5 M n cluster 
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to -50% for a 10 3 M Q clust er; the scatter on N,- „ increases from -10% to 80% for 
the same cluster mass range (ICervino et all 12002 : Cervino & Luridianall2004l) . These 
uncertainties can be reduced by combining large numbers of same-mass, young clus- 
ters to a cumulative equivalent mass «10 5 Mq, since a 'universal' IMF would be fully 
sampled in a cluster of this mass, and measuring the summed values of the combined 
clusters as L(Ha)/M c /=< L(Ha)>/<M f />. 

Stochastic sampling, when combined with photometric and dust reddening un- 
certainties, can also affect age determinations from multi-band photometry, by 'pol- 
luting' samples of young clusters with the addition of old me mbers or loss of young 
ones, and lead to mass und erestimates by up to a factor ~2 (Chand ar et all 1201 Obi : 
Fouesneau & Lancon 11201 d ). Even for massive clusters, i.e., in the absence of stochas- 



tic sampling, observational uncertainties still affect final age determinations at the level 
of a factor ~2 and masses at the 60% leve l. Even achieving thi s level of 'accuracy' 
requires a full set of U,B,V,R,I photometry (Chan dar et alJl20f0aH bh. The uncertainties 
in the ages can be 'mimicked' by models, and be incorporated in the expectations, as 
shown in the left panel of Figure [TJ where the ratio L(Ha)/M c ; is shown averaged over 
two age bins: 1.5-5 Myr and 2-8 Myr. Conversely, uncertainties in the masses are 
carried though, in our treatment, as observational uncertainties. 



3. Results for the Young Cluster Population of M51a 

The diagnostic plot discussed in the previous section should be applied to the low- 
luminosity, low-SFR galaxies where the deficiency in ionizing photon rate has been 
detected, in order to test whether such deficiency may be caused by the 'maximum star 
mass - cluster mass' correlation hypothesis discussed in the Introduction. However, 
reliable identification and, especially, photometric measurements of the low-mass star 
clusters expected to be present in such small galaxies become incresingly difficult be- 
yond «1 Mpc, with ground-based facilties (most N ;o „-deficient galaxies are located 
beyond 1 Mpc). In general, multi-wavelength, at least U-to-I, photometry from the 
HST is required, but such observations are not currently available for those galaxies. 

Age and mass measurements of clusters are now becoming available for more 
massive galaxies, where significant populations of star clusters are present and are 
being investigated to understand the formation and evolution o f these systems (e.g., 
Bast ian et all 2005; Whit more et al. 120101 : IChandar et al.ll2010al fb). These populations 



can be used a test-beds for the validation of the method described in section 2. We 
chose here the cluster population of the well-studied nearby galaxy M51a (8.4 Mpc), 
for which multi-band H ST images are available; this galaxy has super-solar metallicity 
dMoustakas et alj|2010h . an important piece of information, since the ratio L(Ha)/M c / 
depends on metallicity (Figure [T] left panel). 

The right panel of Figure [T] shows the location of the young, <8— 10 Myr, M51a 
star clusters in the L(Ha)/M c ./-versus-M f / plane, where the Ha luminosity has been 
corrected for dust attenuation, and each mass point is the sum of the number of clus- 
ters indicated in the figure. The young clusters are divided into two mass bins (star 
symbols), and the symbols represent: (a) sums made with all clusters whose colors in- 
dicate young ages (light grey symbols), and (b) sums made with only clusters which 
have colors of young populations and also have Ha detections above 3 cr (dark grey 
symbols). 
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Figure 1. LEFT: The expected ratio L(Ha)/M c / as a function of the cluster 
mass M c / both for a universal IMF (crosses a nd asterisks), and for a cl uster-mass- 
dependent IMF (filled squares and tria ngles. IWeidner & Kroupa 2006). The mod- 
els are from Starburst99 (2007 Update, Leithere r et al.ll 19991) . and assume averages 
over large numbers of clusters. The mean L(Ha)/M c / ratios are averaged over two 
age ranges: 1.5-5 Myr (crosses and squares) and 2-8 Myr (asterisks and triangles), 
joined by vertical bars, under the assumption of constant star formation for each age 
range. We assume that clusters younger than 1.5-2 Myr are too dust-enshrouded to 
provide significant contribution at optical wavelengths. The black symbols are for 
1.5 solar metallicity models; examples of higher (~3.5 solar, dark grey) and lower 
(~l/3 solar, light grey) metallicity are shown at a representative mass. RIGHT: The 
ratio <L(Ha)>/<M c i > as a function of the mean cluster mass M c ; in two mass bins 
for M51a, for clusters younger than ~8— 10 Myr (stars with 1 cr error bars), including 
all clusters within each mass bin (light grey, bottom symbols), and including only 
clusters with Ha flux detections >3 a (dark grey, top symbols). The number of clus- 
ters in each mass bin is shown close to the symbols. Expectations from models with 
1 .5 solar metallicity and averaged over the 2-8 Myr age range are shown (asterisks 
and triangles), as they are the closest match to the M51a high-mass data. 
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The data show that in the case of the cluster population of M51a there is no ob- 
vious differential trend in L(Hor)/M c ; for the two mass bins; they are consistent with 
a universal, stochastically-sampled, IMF, although the 'maximum star mass - cluster 
mass' correlation cannot be completely excluded at this stage (at the »2 o~ level). This 
is not the most stringent test of the stochastic sampling model because clusters more 
massive than 10 3 M are expected to contain more than 1 O-type star. 

The data follow the 2-8 Myr age average points (asterisks) of the 1 .5 solar metal- 
licity stellar population models, implying the following, non-exclusive, possibilities: 
(1) our criteria select clusters fairly uniformly distributed in the age range 2-8 Myr; (2) 
the models used in Figured] right, have too low metallicity for M51a, and a higher 
metallicity model should be used (cf. Figure [TJ left panel); (3) significant, mass- 
independent, ionizing photon leakage is present in our regions, which we have not 
attempted to correct for in our analysis. 



4. Conclusions 

While direct measurements of the stellar IMF, even at the upper end where stars are 
bright, become a daunting and ultimately impossible task beyond distances of tens or 
hundreds of kpc from our Galaxy, integral properties of massive stars in coeval popula- 
tions (young star clusters) can be exploited to at least place limits on variations of the 
upper-IMF up to distances of at least 10 Mpc. Specifically, a diagnostic plotting the 
ionizing photon rate per unit star cluster mass as a function of the cluster mass itself 
can provide sufficient discriminating power for some models of IMF variation. 

This method is made more attractive by the increasing availability of well-characterized 
cluster populations in galaxies, thanks to the angular resolution and multi-band cover- 
age of the HST, although there is still dearth of information for the galaxies where a 
deficiency of ionizing photon rate has been observed. 

A number of uncertainties in the age and mass determinations are still present, 
however, from degeneracies in the multi-band colors (which require at least coverage 
in 5 bands from U to I, and possibly UV) and from the stochastic sampling of the IMF 
for cluster masses below «10 4 M©. These will need to be carefully quantified. 

Despite all limitations, some conclusions can be drawn at least for the young clus- 
ter population in the nearby galaxy M51a. The HST data show that the ionizing photon 
rate from these clusters remains constant (when divided by cluster mass) for decreasing 
cluster mass. This suggests that the clusters in M51a are consistent with a univer- 
sal, stochastically-sampled, IMF. This is still a ~2 cr result, and analyses of additional 
galaxies covering a large range of parameters and cluster masses lower than ~10 3 M Q 
will be needed to place this result on firmer ground. 
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